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Abstract 

We show that geometric scaling is satisfied to good accuracy in the forward region of 
d+Au collisions at RHIC. Scaling violations do show up, however, at mid-rapidity, 
and the anomalous dimension of the small-x gluon distribution evolves to near its 
DGLAP limit for transverse momenta of a few GeV. This represents a first con¬ 
sistency check of RHIG deuteron-nucleus and HERA DIS phenomenology, and of 
the universality of the underlying Golor Glass Gondensate (GGC) theory, which de¬ 
scribes both phenomena. It also reconciles successful leading-twist LO and NLO 
perturbative QCD computations of mid-rapidity particle production with small-x 
evolution. Finally, we introduce a new parameterization for the anomalous dimen¬ 
sion of the small-x gluon distribution which properly reproduces known theoretical 
limits at large rapidity, at large virtuality, and on the saturation boundary, and 
still fits the available data from d-|-Au collisions at RHIG. We find indications that 
sub-asymptotic terms in the rapidity-evolution of the anomalous dimension are large. 


1 Introduction 


The weak-coupling gluon saturation formalism nia has been used at HERA, rather suc¬ 
cessfully, in order to describe the x and dependence of the inclusive proton structure 
function F 2 {x,Q‘^). Specihcally, the simple Golec-Biernat and Wiisthoff (GBW) model [Sj 
provided a very efficient qualitative “summary” of the small-x, moderate data in terms 
of an initial condition Q^(xo) = for the saturation momentum and its (constant) growth 
rate A = d\ogQl/d\ogl/x. It led to the discovery of the interesting phenomenon of “geo¬ 
metric scaling” j 3 ], implying that in the relevant range of x and Q^, the structure function 
depends on x and exclusively via the scaling variable Q‘^/Ql{x). In other words, that 
the scattering cross section for a dipole of size r* does not depend on r* and x separately 
but only through the combination p = rtQs{x). 

This phenomenon does not hnd a natural explanation within the linear Dokshitzer, 
Gribov, Lipatov, Altarelli, Parisi (DGLAP) [S] and Balitsky, Fadin, Kuraev, Lipatov 
(BFKL) U QGD evolution equations in the absence of saturation boundary conditions. It 
has been shown to arise from the BFKL equation with saturation boundary conditions [7j , 
in an expansion of the LO-BFKL solution to first order about the saturation line, and 
in mean-held approximation. That equation, with the corresponding boundary condition, 
represents an approximation to the full non-linear “Golor Glass Gondensate” (GGG) theory 
of QGD evolution. It is valid only on one side of the saturation line Qs{x), in the dilute 
regime. 

Moreover, refs. [Zj also computed the scaling violations^ which emerge away from the 
saturation boundary at S> Ql{x). They arise from the diffusion term in the expansion 
of the LO-BFKL solution to second order about the saturation saddle-point, which shifts 
the anomalous dimension 7 of the target gluon distribution by Ay rsj y 1 \og{l/vtQsiy)) 
as one moves away from the saturation boundary {y = \ogl/x denotes rapidity). By 
Q%i.x) r'U (5s(a:)/A^, scaling violations grow to order one (A is a non-perturbative scale in 
the infrared, of order Aqcd); Qgs is the upper boundary of the so-called geometric scaling 
window. Beyond this virtuality, the target gluon distribution is described for example by 
the DLA form, if y and log are sufficiently large. 

A £t to the HERA data in the (“extended”) geometric scaling window Ql{x) < < 

has been performed by lancu, Itakura and Munier (HM) in ref. jS] (see also jUj), 
employing a dipole parameterization which represents an approximate solution to the 
JIMWLK equations cn governing the mean-held evolution of color dipoles. They hnd 
evidence for scaling violations in the HERA data, as predicted by the GGG theory: a 
substantial increase of 7 with l/r^ is required to improve the agreement with the data as 
compared to the original GBW dipole model. 

In recent years, new data on single-inclusive hadron production at semi-hard transverse 
momenta emerged from the Relativistic Heavy-Ion Gollider (RHIG). Here, we are interested 
in particular in data from deuteron-gold collisions at central and forward rapidity (toward 
the fragmentation region of the deuteron projectile) E], which may probe 

^We shall use this terminology throughout the manuscript to refer to violations of geometric scaling 
away from the saturation line, not to the violation of scaling due to DGLAP evolution. 
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high gluon densities in the gold target, and at the same time, should be less affected by 
hnal-state effects than nucleus-nucleus collisions. 

A hrst semi-quantitative analysis of the d-l-Au data from RHIC was performed by 
Kharzeev, Kovchegov and Tuchin (KKT) in ref. ^3], identifying qualitative features of the 
expected CGC evolution and providing an alternative dipole parameterization: the KKT 
model, to be described below, was constructed to £t RHIC instead of HERA data. In 
ref. jin], we extended the standard CGC approach to particle production to include recoil 
effects of the large-a: projectile partons, and found very good quantitative agreement with 
RHIC data at large rapidity using either the KKT or the HM dipole. It turned out that a 
crucial feature of both dipole models lies in the fact that the target anomalous dimension 
at large rapidity and moderate pt is essentially constant, 7 ~ 0 . 6 , and substantially smaller 
than in the DGLAP regime {'Jdglap ~ !)• Hence, for transverse momenta of a few GeV 
and rapidity the data points are in either the saturation or extended scaling regions, 

where scaling violations are essentially absent at RHIC energy. 

In the present paper, we extend our previous calculation of forward hadron produc¬ 
tion to the mid-rapidity region. Our main objective is to analyze whether scaling 
violations, as predicted by the CGC theory and previous dipole hts to HERA data, show 
up. This is a necessary condition for establishing the CGC as a universal theory of QCD 
evolution near a saturation boundary within the weak-coupling (semi-hard) regime. 

The paper is organized as follows. In the next section, we introduce the formula for 
the single-inclusive hadron production cross section in high-energy proton (or deuteron) 
collisions with a heavy nucleus which accounts for both small-x evolution of the dense 
target and full DGLAP evolution of the dilute projectile. We also present the KKT dipole 
parameterization there. In section jS] we apply this formalism to hadron production at mid¬ 
rapidity and compare to data from RHIC to check for the presence of scaling violations. In 
section m we modify the HM dipole parameterization, in particular that of the anomalous 
dimension, to fit the available data from RHIC at central and at forward rapidity, which 
is done in section 0 Finally, section 0 contains our summary and conclusion. 


2 Proton-Nucleus Collisions within the CGC ap¬ 
proach 


The single-inclusive hadron production cross section in p-|-A collision is given by 

d^pA^hX I 


Xf 


/ dxp ^ fg/p{xp, Q)) Nf [^Pu b) Dh/g Q}\ + 

fxp 'i'F XJiF / \'^P / 


dxF d?Ptd?b (27r)2 

DuiA—,Q 

\XF / \ Xy) 


( 1 ) 


where pt and xf are the transverse momentum and the Feynman-x of the produced hadron, 
respectively. Xp denotes the momentum fraction of a projectile parton and b is the impact 
parameter. Note that this expression is different from the more common kt factorized 
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expressions pini in that the radiation vertex is exact. The common eikonal (soft glnon) 
approximation made in the fc^-factorization approach is omitted, which reflects in the fact 
that the projectile parton distribntion fnnctions f{xp,Q^f), and their fragmentation fnnc- 
tions D{z,Q'j) into hadrons, evolve according to the fnll DGLAP evolntion eqnations 
and respect the momentnm snm-rnle. The importance of this evolntion for the interpreta¬ 
tion of forward-rapidity data from RHIC is discnssed in detail in HE]. 

Eq. (P) resnms large logarithms of 1 /xa arising from evolntion of the target wave 
fnnction with energy as well as logarithms of from evolntion of the collinearly factorized 
projectile qnark and glnon distribntions with Q^. It does not inclnde large-a: effects snch 
as recoil in the nncleus, nor does it acconnt for logarithms of on the target side (only 
collinear emissions from the projectile are inclnded, c.f. ref. for details). Therefore, 
this expression is valid when the projectile partons are in either the full or double-log 
DGLAP kinematics (in other words, above the “extended” geometric scaling window for 
the projectile) while the target is in the small-x^, high-density regime. 

To illustrate the regime of applicability of eq. (HD from another perspective we compare 
to the usual fct-factorization approach, where the single-inclusive gluon production cross 
section is given by 

d(j f ri^ 

The pre-factor F = AnasNc/— l)/pt and (pp^A^Xp^A, k^) = dxp^Ad {xp^ a, k^) / {d log k'f) 
denote the unintegrated gluon distribution functions of the proton and of the nucleus, 
respectively. 

To leading logp^ accuracy we only need to pick up the contributions from /c* ~ 0 and 
kt ~ Pt- 

Tdy^t ^ -fo * k^)(j)A{xA,Pt) + * dk^^^p{xp,Pt)d)A{xA, k^) ( 3 ) 

fPt 1 

= Xpg{xp,p^)(j)A{xA,Pt) + M^P^Pt) dkl-^(t)A{xA,kl) . ( 4 ) 

In the second term, we shifted {kt — ptY k'^- 

The hrst term corresponds to eq. (^Q), if one replaces the Dh/g{z) fragmentation function 
from there by 5(1 — z) (we are now discussing gluon, not hadron production) and identihes 


F (l>A{Pt) 


NA{Pt) 

(2vr)2 


( 5 ) 


Nevertheless, one should keep in mind that eq. o is derived in the recoilless approximation 
and so is not applicable when Xp is large, say > 0.01 for a proton projectile. As already 
mentioned above, eq. (HD was derived without this limitation Hi- 

On the other hand, the second term from © is clearly not part of eq. ©. For Pt 
Qgs (the geometric scaling momentum of the nucleus) the unintegrated gluon distribution 
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functions are given by (see e.g. CBI) 


7 




(i)A{xA,kl) ~ SAkj 




( 6 ) 


Here, Sa denotes the transverse area of the nucleus and 7 < 1 is the anomalous dimension 
of its gluon distribution function, see below. Constant prefactors which are irrelevant for 
the present discussion have been dropped. Using these expressions in dSD shows that no 
logarithms of pi arise from the second term, which therefore is suppressed. On the other 
hand, for pt S> Qgs{,XA)i the integral over kl in the second term of (jH) does give rise 
to a contribution log(p 2 /Q 2 j since (j)A{xA,kl) 1 for kt > Qgs- Hence, for pt ^ Qgg 
collinear emissions from the target have to be included as well, for example by using CCFM 
unintegrated gluon distributions m- 

The purpose of the present manuscript is to check for the presence of geometric scaling 
violations in the mid-rapidity data from RHIC. Consequently, we focus on the regime of 
moderately large transverse momenta not too far above Qgg. Since there Xp^ 0.01, we 
believe that it is more appropriate to employ eq. ([Q) rather than ©. 

The cross section depends on the scattering probability of dipoles of size in the 
fundamental and adjoint representations, respectively, at an impact parameter b: 

NfM) = ^Tre(l-U^(6-n/2)U(6 + n/2)), 

NA{n,h) ^ _^Tr, (l-f/t(5_^^/2)f/(6 + r,/2)) , (7) 

where V and U denote Wilson lines along the light cone ^01 in the corresponding repre¬ 
sentation, and Nc is the number of colors. 

In practice, the ab-initio determination of dipole profiles by solution of the JIMWLK 
equations is not feasible yet. Therefore, one resorts to theory-motivated phenomenological 
parametrizations of the dipole profile which can be tested against data. An example is the 
KKT model jT3], where the dipole profile is parameterized as 


NA{rt,yh) 


1 — exp 





( 8 ) 


with Hh the rapidity of the produced hadron. The saturation scale of the nucleus at yh is 
given as^ 


Qs{yh) = Qo exp[A(?//, - yo)/2] . (9) 

Here, A ~ 0.3 is hxed by the DIS data and the initial condition Qo — 1 GeV is set at 
j/o = 0.6, where small-x quantum evolution becomes effective. The saturation scale for the 
dipole in the fundamental representation differs by a factor of Ql —> Ql Cf/Ca = ^Ql- 

^Note that if hadron masses and intrinsic transverse momenta from fragmentation are neglected, then 
the rapidity of the produced hadron equals that of its parent parton. Hence, we do not distinguish between 
them. 
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7 denotes the anomalous dimension of the target gluon distribution with saturation 
boundary condition and is modeled by KKT as 


livh, qt) 


1 

2 



_ \^{yh,qt)\ _ 

\^{yh,qt)\ + ^J2\^{yh,qt) \ + 7cC(3) 


( 10 ) 


Here, we have written 7 as a function of transverse momentum rather than dipole size by 
evaluating it at a characteristic value r* ~ 1/qt- The free parameter c in (jlOj) governs the 
onset of the quantum evolution regime, and has been hxed by KKT to c = 4. 

The Ansatz (jHI) exhibits geometric scaling when 7 is constant. Scaling violations are 
introduced into the KKT parameterization through the function 


^{yh,qt) 


^osjq^m 

{X/2){yh-yo) ' 


( 11 ) 


^ vanishes as yh ^ 00 at hxed qt, so that 7 —1/2. On the other hand, approaching 
the saturation boundary at large rapidity far from the initial condition (i.e. qt Qs{yh) 
for yh ^ yo) leads to ^ 2 . Due to the last term in the denominator of eq. (fTT]|l this 

corresponds to 7 ~ 0.53, again close to the usual BFKL saddle point without saturation 
boundary conditions. 


3 Scaling violations in the mid-rapidity RHIC data 

It was shown in ref. |TH] that the BRAHMS data m on forward hadron production from 
d+Au collisions at RHIC are described very welH by both the KKT ^2] and the HM [S] 
dipoles, provided that DGLAP evolution of the dilute projectile with exact splitting func¬ 
tions is taken into account. The predictions made in appear to be conhrmed by the 
preliminary data from STAR on 7 r° production at rapidity yh — A ^2] • The fact that both 
dipole parameterizations generically include an anomalous dimension of the target gluon 
distribution that differs from 1 at small xa proved crucial. In particular, for yt^ 3 both 
KKT and HM dipoles predict a nearly pt-independent anomalous dimension 7 0.6 for 

Pt^5 GeV, which is close to the LO-BFKL prediction^ of 7 = 7 * ~ 0.627 and, with eq. o 
from above, is consistent with the experimental pt distributions [ini- 

In the present paper, we consider particle production in the central region, yh = 0. As 
we shall see, here the behavior of 7 required by the data and the GGG approach is very 
different. Before comparing to data, however, it is useful to analyze the kinematics. While 
forward hadron production probes rather small xa in the target nucleus (c.f. appendix B 
in ref. cni), one may wonder to what extent the conditions for applicability of the GGG 
formalism are met in the central rapidity region, for RHIG energy [s/s = 200 GeV). We 

^Already at leading order in as, if a phenomenological pt-independent A'-factor for NLO corrections is 
allowed for. 

^This is the anomalous dimension for BFKL evolution with saturation boundary condition, i.e. for 
evolution along the saturation line. 
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expect that the CGC description breaks down beyond some Xq, where recoil effects affect 
the evolution with log 1/x. Near xq^ knowledge of the width Ax/xq of that transition region 
would also be required. Neither xq nor Ax can be determined from the CGC formalism 
itself, since we presently lack a treatment of recoil and of other “large-x” effects within 
this approach. 

However, since evolution in x appears to be quite rapid, it is probably not a bad ap¬ 
proximation to consider a sharp boundary xq beyond which the CGC small-x formalism 
will break down. Experimentally, it is known that HERA data on proton structure func¬ 
tions at or below xq ~ 0.01 can be described reasonably well within the small-x evolution 
approximation Em El 13- In case of nuclear targets, larger values of xq are expected due 
to the ~ enhancement of the saturation scale near xq PI Therefore, for a rough 
estimate, we take Xq = 0.05 as the upper limit of validity of the small-x approximation for 
the target. 

At Uh = 0, the projectile and target x are equal and given by x = qt/\fs^ with qt the 
transverse momentum of the produced parton. If we assume that the small-x approximation 
of the CGC can be applied for x < xq — 0.05, then this translates into transverse momenta 
up to qt — 0.05-\/i. Thus, for RHIC energy the approach should be valid for qtA 10 GeV. 
Hadronization reduces this number by a factor of ~ 0.7, according to the typical values of 
z = xp/Xp m. the fragmentation. Hence, we estimate that in the central region the small-x 
regime extends roughly up to hadron transverse momenta of pt ~ 7 GeV. 



Figure 1: Schematic “phase diagram” of pA collisions at RHIC. For this plot, we have 
simply assumed that Qsiy) = 1 GeVx exp(Ap/2) and Qgsiy) = QK2/)/200 MeV; see, 
however, sections El and El 


This upper limit for the small-x approximation is above the boundary between the 
extended scaling and the dilute DGLAP regimes, Qgs{y) Ql{y)/^ Cl- With Ql ~ 
1 GeV^ at mid-rapidity (averaged over the transverse plane) and A = 200 MeV, this 
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boundary occurs at a parton transverse momentum of about ~ 5 GeV (or hadron 
transverse momentum pt ~ 3.5 GeV). By this qt, the anomalous dimension of the target 
gluon distribution should already have approached its DGLAP limit, starting from 7 = 7 s 
at g* = Qs ~ 1 GeV. Hence, we expect large scaling violations in the central region of p+A 
collisions at RHIG already for transverse momenta on the order of a few GeV, signaling the 
approach towards the DGLAP regime. These different kinematical regimes are summarized 
in Figure lU According to the discussion in the previous section, our approach via eq. ((T)) 
is valid up to g* ~ Qgs but not much beyond. 

These estimates are in qualitative agreement with the results of Accardi and Gyu- 
lassy 121]. They hnd that a resummation of the Glauber multiple-scattering series (which 
basically corresponds to semi-classical saturation) is required at Pi ~ 1 GeV, and that a 
rapid transition to the DGLAP regime occurs towards higher transverse momenta. Gom- 
paring their approach to data from RHIG, they do not seem to hnd evidence for a broad 
geometric scaling window. In this vein, we also recall that a purely semi-classical model 
without quantum evolution does not reproduce the RHIG mid-rapidity data, as shown 
in 1221 ^ (for an alternative approach, see | 23 ])- 

Similarly, leading-twist NLO pQGD m provides an accurate description of the mid¬ 
rapidity cross-section for 5 GeV. We point out that while this indeed contradicts the 
original GBW model, it is perfectly consistent with improved HERA hts including scaling 
violations jH] and with the GGG theory, which does predict that 7 —> ~ 1 for g* > Qgg. 
That this is perhaps the correct interpretation can be inferred from the suppression of 
scaling violations at large rapidity (and pt ~ few GeV), which again is a generic prediction 
of the GGG approach®, and has already been verihed to be consistent with RHIG data [TT)] . 

In Fig. 121 we show the transverse momentum spectra of charged hadrons at mid-rapidity 
for minimum-bias d-|-Au collisions at RHIG^. We also show the theory curves obtained from 
eq. o, using GTEQ5-LO [22] parton distribution functions for the projectile deuteron, and 
KKP-LO j2Zj fragmentation functions at the scale Qf = pt- We plot curves for various fixed 
anomalous dimensions from 7 = 0.5 to 0.9. Glearly, a rather steep rise of the anomalous 
dimension in the pi ~ 1 — 4 GeV range is required for a satisfactory description of the 
mid-rapidity data. At the lower end, pt ^ 1 GeV, the anomalous dimension is close to 
0.6, perhaps marking the onset of saturation dynamics in the weak-coupling regime. This 
observation is further supported by the fact that saturation based models successfully 
describe the hadron multiplicities at RHIG, which are dominated by soft hadrons ^Tj. On 
the other hand, at pt = 4 GeV, the anomalous dimension 7 0.9 required to ht the data 

nearly equals its leading twist DGLAP limit of 1. 

To illustrate this further, we plot the transverse momentum dependence of 7 for the 
KKT model in Fig. (21 For large transverse momenta, we assume that the Fourier transform 
of Va(d) is dominated by transverse distances of order l/g* ~ l/{2pt) and hence evaluate 

fit is possible if one assumes an ad-hoc large “non-perturbative” contribution to Qs, which has to 
be dropped again in the forward region. A critical discussion of this issue is given in ref. |22|. 

^Deriving from the fact that Ay ~ 1/y at fixed pt- See also refs. m 

^For better visibility, we show only the BRAHMS and STAR measurements. Data at central 
rapidity has also been obtained by the PHENIX m collaboration. 




Figure 2 : Charged hadron pt spectrum obtained from eq. (HD using the KKT dipole for 
various hxed anomalous dimensions 7 . The points with error bars show the BRAHMS HU 
and STAR data for minimum-bias d-l-Au collisions (RHIC, mid-rapidity). 


7 (rt) at rt = l/( 2 pt). 

One observes that the KKT parameterization of the anomalous dimension does not vary 
much with transverse momentum, staying within 20 % of 7 ~ 0 . 6 . As discussed above, 
this behavior is generic to the CGC at large rapidity (i.e. small xa), and is consistent with 
the RHIC data Hi- However, we have also seen that the data requires a rapid approach 
towards the DGLAP regime at mid-rapidity and ~ 4 GeV, while 'Jkkt 1 only at 
asymptotically high pt. In other words, we hnd evidence for substantial scaling violations 
at mid-rapidity, larger than those incorporated into the KKT parameterization but in 
agreement with our estimates on Qgg from above, which then weaken at large rapidity 
(they get pushed to higher pt). Hence, the qualitative behavior of 7 implied by the RHIC 
data is remarkably consistent both with HERA phenomenology and with the predictions 
of the CGC approach. 

4 New parameterization of the anomalous dimension 

Motivated by these results, we introduce a new parameterization of the dipole prohle with 
an anomalous dimension incorporating stronger scaling violations than the KKT dipole. 
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Figure 3: Anomalous dimension 'y{yh,rt) of the small-x nuclear gluon distribution for the 
KKT dipole model, evaluated at r* = l/{2pt). 


At the same time, it corrects some dehciencies of the IIM dipole at large Q^. In our 
parameterization, the anomalous dimension is given by 


Ay 


7 , + A-f{Q^,Y) 

(1 _ ^ ) log(QVQf) 

AF + log(QVQ2) + ^^ • 


( 12 ) 


We have written 7 in terms of the variables appropriate for DIS, i.e. and Y = log(l/x), 
similar to the GBW and IIM models. Here, = 1/r^ is the inverse transverse size of the 
dipole, not to be confused with the factorization scale which enters the parton distribution 
and fragmentation functions in eq. ([Q). We shall transform into the variables appropriate 
for hadron-hadron collisions shortly. 

The function is constructed such as to satisfy the following limits: 

1 . at any hxed rapidity H, 7 ^ 1 for —»• 00 . 

2 . if F —>• cx) at hxed Q'^/Q^iY), Ay decreases asymptotically like ~ 1 /F, as determined 
by BFKL jS]. This indicates that the geometric scaling window broadens with Y. 

3. 7 —»• 7 s for —>• Ql{Y) at any rapidity. 

4. At large but hxed rapidity, Y ^ (d/A)^, the geometric scaling window reaches up to 

\ogQ%{Y)/Ql{Y) rsj XY, consistent with the asymptotic estimate for the geometric 
scaling window 0 : ~ Qt/AlcD- 


10 





























This last condition is the most essential difference between our new parameterization and 
the IIM parameterization [S], where 


b') 


7 .+ 


log(QVQg) 

nXY 


(13) 


Here, one hnds \ogQ'^g/Ql (1 — 'ys)K,XYi with k ~ 10. Hence, at large Y, their Q'^g{Y) is 
much larger than Q^^iY)/h?. Such behavior arises if the geometric scaling window is esti¬ 
mated from the diffusion term in the expansion of the LO-BFKL solution to second order 
about the saturation saddle-point. On the other hand, if Qg^ is estimated via the transition 
point between the LLA and DLA saddle points, respectively, one hnds logQ'^^/Q^ ~ XY, 
without the additional factor of k (for more details see, for example, section 2.4.3 in ref. [2] 
and references therein). 

With kXY replaced by XY, however. Ay rises too rapidly with logQ^ to allow for a good 
description of the RHIC data. Also, one should keep in mind that the IIM parameterization 
in fact provides a good ht to the HERA data on F 2 {x,Q^) at small x. Hence, strong 
modihcations of 7 are not desirable. As a consequence, to restore (approximate) agreement 
with the IIM parameterization, we are led to introduce terms which are subleading in Y. 
This agrees with indications from studies of small-x evolution equations which show that 
at realistic energies (rapidities), subleading corrections to the asymptotic expressions are 
quite large. In (HI, these are modeled by the a/R term with a free coefficient. From the 
mid-rapidity RHIC data, we extract d ~ 1.2 (see next section). 



Figure 4: Anomalous dimension versus Q at various x. 

Fig. m compares our dipole parameterization m with d = 1.2 to that from eq. m- 
At small X they are very similar as 7 becomes flatter and approaches 7 <j. Towards large x, 
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7 is somewhat smaller for our parameterization®. At the same time, the new anomalous 
dimension is closer to 1 than 'yxKTi which was shown in Fig. |21®, i.e. our parameterization 
features a narrower scaling window. With d = 1.2, which restores agreement with the IIM 
dipole at a; = 10“® and also describes the RHIC d+Au data correctly, the subleading term 
in the rapidity evolution of 7 actually dominates until Y ~ (d/A)^ = 16. 

Finally, we also point to another recent parameterization of the dipole prohle in ref. p, 
which is meant to reproduce the RHIC data. That parameterization also includes breaking 
of geometric scaling at mid-rapidity, but not the sub-asymptotic terms for the evolution 
with rapidity. The dipole prohle at large rapidity again resembles the KKT, IIM, and 
our present dipole model and so the forward hadron spectra “generated” by all these 
parameterizations should be rather similar. 


5 Application to d+Au collisions at RHIC 

For hadron-hadron collisions, we need to replace l/r^ and Fourier transform the 

dipole prohle function from (jS)): 

dVt NA,F{rt) = dn n Joinqt) NA^ri^t) • (14) 

Unfortunately, it is very challenging to perform this Fourier transform numerically when 
7 is not constant but rj-dependent. It is therefore common to replace 7 (rt) by a constant 
7 (ri = 1/qt), which makes the Fourier transform tractable |151 ITT)]. This is a reasonable 
approximation in the region where 7 (rt) is rather hat, e.g. at large rapidity and pt -C Qgs- 
On the other hand, at mid-rapidity, where scaling violations are larger, we found that it is 
important to acconnt for the r^-dependence of the anomalous dimension, as this strongly 
ahects N{qt). 

However, with a r^-dependent anomalous dimension, we were unable to perform the 
Fourier transform m numerically at very large transverse momentum, where the phase 
factor oscillates extremely rapidly. More importantly, the Fourier transform turns negative 
already at intermediate (~ 8 GeV for gluons at mid-rapidity), where our numerical 
methods are still reliable. This unphysical behavior occurs for all of the above-mentioned 
'y{rt) parameterizations (KKT, IIM and ours) and presumably indicates that the behavior 
of the function at very small r* is inadequate We are nevertheless able to determine the 
transverse momentum spectrum of hadrons reliably up to about p* = 4 GeV, which doesn’t 
receive contributions from the unphysical oscillations of NA{qt) at high parton momentum. 
As discussed in sec. El we do not expect that eq. (CD is valid much beyond this pt at RHIG 
(mid-rapidity) anyways. 

®Ref. 13 considered only the region x < 0.01 and the fact that 7 //m is not bounded from above was 
irrelevant for their analysis. 

®Note also the factor of 2 difference in the pt-scale at which 7 is evaluated. 

^'^This problem was encountered also in earlier studies m of Fourier transforms of the Bartels et al. 
DGLAP-improved dipole profile (3. 
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The field of the target nucleus is probed at the rapidity i/a given by ua = logl/a:A = 
log(l/a:p) + 2yh (see appendix B in [12] for details of kinematic definitions). Hence, the 
saturation scale is defined at that rapidity, 

Ql {xa) = Ql A]!^ {xq/xa)^ . (15) 

This implies that at fixed hadron rapidity yh, as we sum the contributions from projectile 
partons with different momentum fractions Xp in eq. o, Qs{xa) also varies. 

The effective mass number of the nucleus depends on the impact parameter; following 
ref. pnj, for heavy A ~ 200 targets we do not perform the integral over b explicitly but 
take Hgff = 18.5 for minimum bias collisions (the issue of impact parameter averaging is 
also discussed in ref. |22|)- We fix the scale xq = 3 ■ 10“"^ to the value extracted in [Sj 
by a fit to HERA data. With Qo = 1 GeV, one obtains Qs — 1 GeV at mid-rapidity for 
transverse momenta ~ 1 GeV, which is reasonable ng. 



Figure 5: Gomparison of theory and data uniiaiii for minimum-bias d-|-Au collisions at 
RHIG energy. 

In Fig. 121 we show our results for minimum-bias d-|-Au collisions at RHIG energy with 
the new parameterization m- The factorization scale in eq. o is taken as Qf = pt for 
all rapidities. All other parameters have been fixed above. As promised, the shape of the 
mid-rapidity pf-distribution is described quite well. For the correct overall normalization, 
we need to multiply our LO calculation by a factor Ky=Q ~ 3. This is much larger than at 
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forward rapidity, were ~ 2 and Ky= 4 ^ ~ 1 , respectively (for the same scale Qf = Pt 

and the same parton distribution and fragmentation (27j functions). Such a rapid 
increase of the i^-factor towards mid-rapidity is expected because of larger phase space for 
NLO corrections Despite the good agreement of the LO curve with the prdependence 
of the data, NLO calculations of particle production are clearly desirable. 

Fig.ISlalso shows that the new parameterization of the anomalous dimension introduced 
here describes the forward data as well as the IIM and KKT dipole models (for the latter, 
see curves in uni). As already alluded to above, the reason for this generic behavior is that 
all models predict a rather flat 7 (rt) 7 ^ at small x (wide geometric scaling window), and 

hence do not differ much. This also reflects in the fact that the forward rapidity curves 
from Fig. El which were obtained with the r^-dependent anomalous dimension, are quite 
similar to those shown in ref. which employed a constant 7 = 7 ( 7 "^ = l/2pt). We 
stress that such good agreement with the data in the forward region is only achieved if the 
projectile parton distributions satisfy DGLAP evolution with exact rather than small-x 
approximated splitting functions |lfij . 

6 Summary 

In summary, in this paper we extended our previous studies of forward rapidity hadron 
production in deuteron-gold collisions at RHIC to the central rapidity region. We argued 
that here the CGC predicts signihcantly larger scaling violations than at forward rapidity, 
and that indeed a relatively rapid transition of the anomalous dimension of the target 
gluon distribution from 7 ^ 0.63 (BFKL with saturation boundary conditions) to 7 —> 1 

(DGLAP) is required by the data. It is important to note that this transition occurs in a 
regime where the small-x approximation made in the Golor Glass Gondensate formalism 
is still valid. Hence, the scaling violations predicted by the GGG can indeed be tested. 

We propose a new parameterization of the dipole prohle which features a steeper trans¬ 
verse momentum dependence than the KKT dipole and leads to a satisfactory description 
of the mid-rapidity data. This agreement supports the predictions of the Golor Glass 
Gondensate formalism regarding the existence of different kinematic regions with different 
underlying physics. 

Schematically, at small transverse momentum one probes the saturation region of the 
nucleus, where the underlying physics is that of a “black” target. As the transverse mo¬ 
mentum increases, there is a transition from the saturation to the (extended) geometric 
scaling regime, corresponding to an anomalous dimension close to that from BFKL with 
saturation boundary condition. A further increase in the transverse momentum takes one 
beyond the geometric scaling region, signaled by the necessity of introducing rather large 
scaling breaking in the dipole prohle. We showed that for transverse momenta of about 
4 GeV, one is already in a small-a; DGLAP regime. Thus, the transition from the satura¬ 
tion to the DGLAP regions is quite narrow in the central rapidity region at RHIG. On the 
other hand, at large rapidity (?/h~ 3) only the saturation and scaling regions remain, while 

^^This trend is also present in pQCD calculations of hadron spectra in proton-proton collisions ESI 
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the DGLAP window is cut off by finite-energy constraints, which has not been realized 
before. At LHC energy, all these kinematic regimes extend to even higher pt, providing a 
great opportunity to study the various “phases” of high-energy QCD in more detail and 
with better accuracy. 
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